Hydrodynamical analysis of centrality dependence of charged particle's multiplicity in 

v^^jv=2.76 TeV Pb+Pb collisions 

A. K. ChaudhurE and Victor Rojill 

Variable Energy Cyclotron Centre, 1/AF, Bidhan Nagar, Kolkata 700 064, India 

(Dated: January 13, 2013) 

In Israel-Stewart's theory of dissipative hydrodynamics, we have analyzed the recent ALICE data 
for the centrality dependence of charged particle multiplicity, centrality dependence of integrated 
elliptic flow and 0-5% charged particles pt spectra in -^5^^^=2.76 TeV Pb+Pb collisions and deter- 
mined the initial or the thermalisation time. Analysis indicate that the ALICE data disfavor very 
early thermalsiation ri=0.2 fm, or very late thermalsation Ti= 4fm. Data are best explained for 



thermalisation time Ti — 0.88 



+0.68 



fm. 
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PACS numbers: 47.75. +f, 25.75.-q, 25.75.Ld 

Lattice simulations of QCD indicate that the strongly 
interacting nuclear matter can undergo a confinement- 
deconfinement cross-over transition Experiments 
at Relativistic Heavy Ion Collider (RHIC), produced con- 
vincing evidences that in y/sj^j,f=200 GeV Au-f Au colli- 
sions a collective QCD medium is created [B-Q, though it 
is uncertain vi^hether or not the matter produced can be 
characterised as Quark-Gluon Plasma (QGP), the lattice 
QCD predicted deconfined phase. The issue is expected 
to be settled at Large Hadron Collider, where Lead nu- 
clei will collide head on at enormous energy ^/s]^J^=5.5 
TeV. Recently, ALICE collaboration measured the 

centrality dependence of charged particle multiplicity in 
yijVAr=2-76 TeV Pb-hPb collision. The centrality de- 
pendence is similar to that obtained in y^^^=200 GeV 
Au-f Au collisions at RHIC, charged particle multiplicity, 
normalized by the participant nucleon pair increase by 
a factor of 2 from peripheral (70-80%) to central (0-5%) 
collisions. ALICE collaboration also published charged 
particles px spectra and elliptic flow [HI, In 0-5% 
collisions, px spectra of charged particle's are suppressed 
by a factor Raa ~ 0.14 at pT=6-7 GeV, which is smaller 
than at lower energies. In peripheral collision, suppres- 
sion is modest, Raa ~0.6-0.7. Integrated elliptic flow is 
~30% more than that in Au+Au collisions at RHIC. 

Hydrodynamic models have been used extensively 
to analyze the experimental data in ^^^^=200 GeV 
Au+Au collisions at RHIC and obtain information about 
the initial condition of the produced medium. Hydrody- 
namic models require the assumption of 'local' thermal 
equilibrium. If the assumption is met, relativistic hy- 
drodynamic equations can be solved to trace-back to the 
initial fluid condition from experimental data. At RHIC 
energy, near ideal QGP fluid, initialized to central en- 
ergy density Si « 30 GeV/ fm^ at initial time n ssO.6 
fm, explains a large variety of experimental data [l3l - [T7| . 

One of the important issue in a hydrodynamical model 
analysis is the initial or the thermalisation time. In the 



present letter, we have analysed the ALICE data for 
the centrality dependence of charged particle multiplic- 
ity , centrality dependence of integrated elliptic flow 
13] and charged particles px spectra in 0-5% collision 
ll|, to determine the thermalisation time in LHC en- 
ergy collisions. The three data sets are best explained 
if thermalisation time is sal fm. Small thermalisation 
time Ti=0.2 fm, or large thermalisation time Ti—A fm are 
disfavored by the data. 

We assume that in \/s^jy=2-76 TeV Pb+Pb collisions, 
a baryon free QGP fluid is formed. Viscosity to ent ropy 
ratio of the fluid is assumed to be ry/s = l/47r [l^, [T9| . 
The space-time evolution of the fluid is obtained by solv- 
ing Israel-Stewart's 2nd order theory [20l. [2]|. 
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EqU is the conservation equation for 
momentum tensor, T^'^ — {e + p)u^u^ — pg'"" + tt'' 
£, p and u being the energy density, pressure and fluid 
velocity respectively, tt^" is the shear stress tensor. Eql2] 
is the relaxation equation for the shear stress tensor ■n^'^ . 
In Eql^l D ~ w^S^ is the convective time derivative, 
y<t^u''> = KV^m" + Vuf") - ^id.u)igt"' - m^m") is a 
symmetric traceless tensor, rj is the shear viscosity and r^r 
is the relaxation time. It may be mentioned that in a con- 
formally symmetric fluid relaxation equation can contain 
additional terms [l^. Assuming boost- invariance, Eqsd] 
and [2] are solved in (r = \/t'^ — z'^,x, y, rjs = ^Inf^f) 
coordinates, with the code " 'AZHYDRO-KOLKATA'' ', 
developed at the Cyclotron Centre, Kolkata. Details of 
the code can be found in 
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Hydrodynamic equations are closed with an equation 
of state p = p{s). In the present study, we use an equa- 
tion of state where the Wuppertal-Budapest 3] lattice 
simulations for the deconfined phase is smoothly joined 
a,t T — Tc = 174 MeV, with hadronic resonance gas EoS 
comprising all the resonances below mass mres=2.5 GeV. 
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FIG. 1: (color online) The black circles are ALICE data for 
centrality dependence of charged particle's multiplicity per 
participant pair in 2.76 GeV Pb+Pb collisions [lol |. 

The black, red, green,yellow and blue lines are hydrodynamic 
simulations with initial time Ti— 0.2, 0.6, 1.0, 2.0 and 4.0 fm 
respectively. 
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FIG. 2: (color online) Black, red, green, yellow and blue lines 
are hydrodynamical model simulations for elliptic flow for pi- 
ous, for initial time Ti=0.2, 0.6, 1.0, 2.0 and 4.0 fm. The black 
circles are the ALICE measurement '12] for charged particle's 
elliptic flow in 4-particle cumulant method. 



Solution of partial differential equations fEas lll2p re- 
quires initial conditions, e.g. transverse profile of the 
energy density {e{x,y)), fluid velocity {vx{x,y),Vy{x,y)) 
and shear stress tensor {n'^'^ {x , y)) at the initial time r^. 
One also need to specify the viscosity (77) and the relax- 
ation time (t^). a freeze-out prescription is also needed 
to convert the information about fluid energy density 
and velocity to particle spectra and compare with ex- 
periment. We assume that in an impact parameter b 
collision, at the initial time r^, initial energy density is 
distributed as in a Glauber model llSll. 



^(b, X, y) = e4(l-/)iVpart(b, x, y) + fNcoii{h, x, y)], (3) 

Npart and Ncoii in Eq|3] are the transverse profile of 
the average number of participants and average number 
collisions respectively. / is the hard scattering fraction. 
We assume that -^5^^=2.76 TcV Pb+Pb collisions are 
dominated by hard scattering and fix / = 0.9. We also 
assume that at the initial time fiuid velocity is zero. We 
initialize the shear stress tensor to boost-invariant value, 
^xx _ ^yy _ 2ri/3Ti, TT^y=0. FoT the relaxation time, we 
use the Boltzmann estimate r^r — Sri/2p. Fluid viscosity 
is assumed to be 77/s = l/47r and we neglect any temper- 
ature dependence of rj/s. Hydrodynamics also require a 
freeze-out condition. We assume that the fluid freeze-out 
at a fixed temperature Tf=130 MeV. In viscous hydro- 
dynamics, particle production has contributions from the 
non-equilibrium part of the distribution function. We in- 
clude the non-equilibrium contribution. We note that 
the freeze-out condition does not account for the validity 
condition for viscous hydrodynamics, i.e. relaxation time 
for dissipative fluxes are much greater than the inverse 
of the expansion rate, trO^u^^ << 1. Recently, Dusling 
and Teaney [23| implemented dynamical freeze-out con- 
dition. In dynamical freeze-out, non-equilibrium effects 



are stronger than obtained at fixed temperature freeze- 
out. Freeze-out at Ti?=130 MeV however satisfies the vis- 
cous hydrodynamic condition that non-equilibrium con- 
tribution to particle production is much smaller than the 
equilibrium contribution. 

The central energy density £i at the initial time is 
fixed to reproduce experimental charged particle's mul- 
tiphcity = 1601 ± 60 in 0-5% Pb-^Pb collision. For 
initial time Ti—Q.2, 0.6, 1.0, 2.0 and 4.0 fm, we have simu- 
lated 0-5% Pb-f Pb collisions and computed negative pion 
multiplicity. Resonance production is included. Noting 
that pion's constitute w85% of the total charged parti- 
cles, TT" multiplicity is multiplied by the factor 2 x 1.15 to 
compare with experimental charged particle multiplicity. 
Irrespective of the initial time, central energy density can 
be varied to reproduce experimental multiplicity in 0-5% 
collision centrality. In tableHl central energy density and 
temperature required to reproduce experimental charged 
particle multiplicity in 0-5% collisions are listed. Uncer- 
tainty in Ei or Ti refiect the uncertainty in ALICE mea- 
surements. If QGP fluid is thermalised in the time scale 
ri=0.2 fm, very large energy density ^ 370 GeV/ fm^ is 
required to reproduce the experimental multiplicity. En- 
ergy density is reduced if thermalisation time increased. 
One also note that for comparable thermalisation time, 
at LHC, fluid is produced at much higher energy density 
or temperature than at RHIC energy collisions. For ex- 
ample, for Ti=0.6 fm, Si « 126 GeV/ fm? at LHC energy 
collisions, as compared to et ~ 30 GeV/ fm? at RHIC 
collisions. 

With the initial conditions fixed, we have simulated 
Pb-|-Pb collisions over a wide range of collision central- 
ities. In Fig [T] simulation results are compared with the 
ALICE measurements for the centrality dependence of 
charged particle multiplicity per participant nucleon pair 
(-5^^^^^). ALICE measurements are not reproduced 
if the fiuid thermalises early, ri=0.2 fm, multiplicity in 
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TABLE I: Central energy density (e^) and temperature (Ti) required to reproduce experimental charged particle's multiplicity 
in 0-5% Pb+Pb collisions, for ri=0. 2-4.0 fm. Freeze-out temperature Tf = 130 MeV. Model estimates for the charged particle 
multiplicity is also noted. In the last 4 columns, /N for the data sets analysed are noted. 
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"data for 0-5% PT spectra are read from [llj . We assume 7% error 
in the data. 



peripheral collisions are over predicted. Data are bet- 
ter reproduced for >0.2 fm. Very different trend of 
the simulated — T^^^ '^«=0.2 fm is puzzling. In 
a Glauber model of initialisation, one expects that with 
different initial times, centrality dependence will be qual- 
itatively similar, though differ in quantitative magnitude. 
Indeed, simulation results do indicate that for initial time 
Ti >0.2 fm, though quantitatively different, qualitatively, 
centrality dependence of multiplicity per participant pair 
is similar. For Ti— 0.2, qualitative nature is changed. 
The reason can be understood as follows: For Ti— 0.2 
fm, initial pressure, even in a peripheral collision, is very 
high. Due to large pressure gradient fluid can accelerates 
rapidly and achieve large velocity. We have checked that 
for Ti=0.2 fm, in central and peripheral collisions, at the 
freeze-out fluid velocity (vr) can be as high as Vr=0.8c. In 
contrast, for larger thermalisation time Ti=.6 fm, though 
in a central collisions, freeze-out velocity is Vr ~0.8c, in a 
peripheral collision, Vr ^0.6c. Even though initialisation 
with Ti=0.2 and 0.6 fm, both produces similar number of 
particles in a central collisions, in peripheral collisions, 
initialisation with Ti— 0.2 fm produces more particles due 
to increased fluid velocity. We may note here that at 
very early time, e.g. Ti=0.2 fm, non-equilibrium effect 
can be large and even viscous hydrodynamics may not 
be applicable. Failure to explain the data with initial 
time Ti=0.2 fm may be either due to improper choice of 
initial time, or due to inapplicability of viscous hydrody- 
namics or both. But charged particles multiplicity data 
definitely disfavor very early initialisation. 

Integrated elliptic flow is an important observable in 
relativistic energy collisions. It is a measure of collec- 
tivity in the medium. In Figj^l centrality dependence of 
simulated (integrated) elliptic flow in -y/sjY^=2.76 TeV 
Pb-t-Pb collisions are compared with ALICE measure- 
ments [111- Even though Ti— 0.2 fm is not consistent 
with centrality dependence of charged particle multiplic- 
ity, simulated flows agree with the experiment. Indeed, 
for ri=0. 2-1.0 fm, simulated flows are nearly identical 
(differ by less than 3%) and well explains the ALICE 
data. For higher thermalisation time > 2 fm, descrip- 



tion to the data gets poorer. Also V2 in peripheral colli- 
sions continue to increase as oppose to the experimental 
data. Naively, one expects that for late initialsiation, 
pressure will be low and expansion time will shorten and 
elliptic flow will not grow. Explicit simulation indicate 
that for late thermalisation, though the fluid pressure is 
low, expansion time is not shortened likewise. For ex- 
ample, for initial time t,;=0.6 fm, in a peripheral colli- 
sion, (central) fluid freezes out at Tp «6.6 fm and for 
Ti=4.0 fm, expansion time is shortened marginally, fluid 
freeze-out at, Tp «9.9 fm. However, growth of momen- 
tum anisotropy, Sp = j darfy(T^=^+T»«) ^i^'^^ markedly. 
For Ti=0.6 fm, Sp, after few fm of evolution saturate to 
Ep ^0.13. Due to decreased pressure, for late initiali- 
sation, momentum anisotropy grow rather slowly. It do 
not saturate, but continue to grow till freeze-out and at 
freeze-out Sp ~0.17. Enhanced momentum anisotropy 
for late thermalisation is reflected as enhanced V2. 

In Figl3]ALICE measurements [llj for the charged par- 
ticle's spectra in 0-5% collision are compared with hydro- 
dynamical simulations. Slope of the spectra is enhanced 
for smaller initial time. It is expected also. Other con- 
ditions remaining unchanged, fluid initialized at higher 
temperature produces more high pT particles than fluid 
initialized at lower temperature. 0-5% spectra are rea- 
sonably well explained in simulations with thermalisation 
time rj=0. 6-2.0 fm. Description to the data gets poorer 
for thermalisation time > 2 fm. 

For a quantitative analysis, we have computed /N 
values for the data sets analysed. They are noted in ta- 
bleHl From the values, it is apparent that if only a 
single data set is considered, one may largely underesti- 
mate or overestimate the thermalisation time. For exam- 
ple, if one consider only the elliptic flow data, one may 
conclude that the thermalisation time at LHC is Ti=0.2 
fm. However, thermalisation time Ti=0.2 fm is not pref- 
ered by the ALICE data on the centrality dependence of 
the charged particles multiplicity. Perfect fit to the mul- 
tiplicity data is obtained for thermalisation time ri=1.0- 
2.0 fm. 0-5% pt spetcra in 0-5% collisions on the other- 
hand prefer thermalisation time Ti=1.0 fm. Thermalisa- 
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FIG. 3: (color online) black circles are charged particle's spec- 
tra in 0-5% Pb-fPb collisions at ^^^^=2.76 TeV The 
black, red, green, yellow and blue lines are hydrodynamic 
model simulations for charged particle's spectra for initial 
time Ti=0.2, 0.6, 1.0, 2.0 and 4.0 fm respectively. 



tion time at LHC energy should be obtained by analysing 
all the data sets simultaneously. In the last column of ta- 
bleHl we have noted the x^/^ values for the three data 
sets combined. The combined data are best expalined 
for thermalisation time Ti—1 fm, (x^/A^)min=2.8. Very 
small thermalisation time Ti— 0.2 fm or very large ther- 
malisation time Ti=4.0 fm is not favored by the data, 
X^/N increase by a factor of ^4-5 from the minimum 
/N value. Using a parabolic fit to the ~x? 1^ values. 



we could estimate the thermsalisation time rather accu- 
rately, Ti = 0.88^g gg fm. It must be mentioned that the 
present estimate is obtained with a specific set of initial 
conditions, e.g. initial zero fluid velocity, hard scattering 
fraction f=0.9, boost invariant values for the initial shear 
stress tensors etc. We have also assumed a fixed viscosity 
to entropy ratio 77/5=0.08. All the possible initial con- 
ditions are not explored. If uncertinty over the initial 
conditions are included, the present estimate of thermal- 
isation time at LHC energy ri=0.74-1.56 fm, will be even 
more uncertain. 

In the present analysis, we have neglected bulk viscos- 
ity. In general bulk viscosity is much smaller than shear 
viscosity. In QCD, near the transition point, bulk vis- 
cosity can be large (H, [2^ . Effect of bulk viscosity on 
particle spectra and elliptic flow has been studied in . 
Compared to shear viscosity, effect of bulk viscosity is 
small. Present result that ALICE data disfavor small or 
large thermalisation time, or the present estimate of the 



thermalisation time t, = 0.88 



-0.14 

+0.68 



fm will also remain 



largely unaltered even if bulk viscous effects are included. 

To conclude, analysing the recent ALICE data for the 
centrality dependence of charged particle multiplicity, in- 
tegrated elliptic flow and pt spectra, we have determined 
the thermalisation time for QGP fluid in 
TeV Pb+Pb collisions. The combined data are best ex- 
plained for thermalisation time = 0.88]^Q gg fm. Small 
thermalisation time rj=0.2 fm, or large thermalisation 
time Ti= 4 fm, is not favored by the data. 
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